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SUMMARY 


OBJECTIVE 

LOWTRAN,  a  computer  code  that  predicts  infrared  properties  of  the  atmosphere,  is 
primarily  intended  for  aerial  and  terrestrial  applications.  LOWTRAN  has  recently  been  modiT  -d 
to  improve  its  predictive  ability  in  marine  eoviromnents.  The  objective  of  this  report  was  to 
describe  a  further  improvement  in  the  marine  modifications.  Sun  glint,  the  specular  reflection  of 
solar  rays  by  the  ocean  surface  into  the  observer’s  line  of  sight,  has  been  included.  The  glint 
model  is  based  on  the  ocean  wave  slope  statistics  of  Cox  and  Munk  (1954). 

RESULTS 

TTie  results  of  the  modified  model  have  been  compared  with  infrared  data  taken  at  3  to 
5  microtrs  over  the  open  ocean  under  various  conditions  of  solar  elevation  and  wind  speed.  A 
significant  glint  contribution  was  obtained  under  the  prevailing  conditions  for  one  data  set: 

1450  PDT  on  10  October  1991.  A  close  match  existed  between  the  shape  of  the  predicted  and 
observed  radiance  values  in  the  ocean  sky  and  sea,  but  the  absolute  value  of  the  predicted  and 
observed  radiances  differed  by  a  large  constant  whose  equivalent  blackbody  temperanire  was  on 
the  order  of  5“C.  The  origin  of  this  additional  radiance  is  unknown.  It  may  be  multiple  scatter¬ 
ing  of  solar  radiation. 
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INTRODUCTION 


In  the  visible  wavelength  range,  nighttime  visibility  is  poor  because  the  primary  source  of 
radiation,  the  sun,  is  hidden.  However,  because  ordinary  objects  on  the  surface  of  the  earth  also 
radiate  energy  both  day  and  night  (albeit  in  the  near*  and  mid-infrared  ranges),  it  is  possible  to 
see  at  night  though  the  eyes  of  an  infrared  camera.  The  spectral  radiance  (mW  cm'^  str*^  lun'^) 
of  a  300-K  blackbody  is  shown  in  figure  1,  where  it  is  compared  with  the  spectral  radiance  of  a 
6000-K  blackbody,  which  is  a  fair  approximation  to  the  solar  spectral  radiance. 

It  is  an  ironic  fact  that  the  departure  of  the  earth’s  atmosphere  from  the  300-K  blackbody  law 
makes  infrared  vision  possible.  This  is  because  the  atmosphere  of  the  earth  (as  well  as  the 
surface  of  the  earth  and  the  surface  of  the  ocean)  has  radiances  that  approach  the  radiance  of  a 
300-K  blackbody,  but  only  at  those  wavelengths  where  the  radiances  absorb  strongly.  Both  the 
earth  and  ocean  absorb  strongly  in  the  near-infrared,  so  their  radiances  are  well  approximated  at 
their  respective  temperatures  near  300  K  by  the  blackbody  function.  The  atmosphere  contains 
two  regions  from  approximately  3  to  5  microns  and  again  from  approximately  8  to  14  microns 
where  absorption  is  low,  especially  for  paths  with  small  optical  depths;  for  example,  absorption 
is  low  direcUy  up  toward  the  zenith  or  horizontally  from  a  high  mountain  range  to  deep  space 
(Bell  et  al.,  1960).  In  these  spectral  regions,  the  radiance  of  the  atmosphere  drops  below  the 
blackbody  value.  This  is  fortunate,  however,  because  then  the  earth  or  sea  (or  a  ship)  can  be  seen 
through  the  relatively  transparent  atmosphere.  Outside  these  regions,  infrared  eyes  are  blinded  in 
the  imediate  vicinity  of  the  camera  by  the  blackbody  radiance  of  the  atmosphere.  The  transpar¬ 
ency  of  the  atmosphere  in  these  spectral  regions  in  the  mid-  and  far-infrared,  together  with 
recent  advances  in  infrared  detector  technology,  make  knowledge  of  the  infrared  radiance  and 
transmission  in  the  vicinity  of  the  earth  important  to  the  armed  services.  This  knowledge  affects 
the  services’  ability  to  operate  at  night  almost  as  well  as  they  can  operate  during  the  day. 

This  document  concerns  radiance  models  that  can  be  used  to  predict  the  infrared  behavior  of 
the  sea,  land,  and  sky.  Personnel  can  compare  the  predicted  behavior  with  experiments. 

MARINE  MODIFICATION  OF  LOWTRAN  6 

The  Atmospheric  Transmittance/Radiance  Computer  Code  LOWTRAN  6  (Kneizys,  et  al., 
1983)  has  recently  been  modified  by  Wollenweber  (1988),  who  introduced  radiance  reflections 
and  blackbody  emissions  from  the  sea  surface  into  the  code  sky  in  order  to  provide  a  more 
realistic  representation  of  the  infrared  marine  environment.  In  this  report,  we  describe  the 
inclusion  of  an  additional  contribution  that  is  also  important  in  marine  settings:  sun  glint,  which 
is  the  specular  reflection  of  solar  rays  into  an  infrared  camera  by  the  ocean  surface. 

The  general  idea  of  the  maritime  modification  of  the  LOWTRAN  6  computer  code  is 
presented  in  figure  2.  The  camera,  located  a  short  distance  above  the  ocean  surface  typical  of 
shipboard  heights  (e.g.,  30  m),  is  adjusted  to  view  the  sky  and  sea  in  the  vicinity  of  the  horizon, 
so  its  zenith  angle  usually  falls  within  a  degree  or  so  of  %  degrees.  The  horizontal  range  will  be 
on  the  order  of  a  few  tens  of  km  at  an  altitude  of  30  m.  Since  typical  instantaneous  fields  of  view 
for  infrared  cameras  are  on  the  order  of  1  mrad,  the  linear  dimensions  of  a  single  camera  pixel  at 
the  ocean  surface — the  pixel  footprint — will  be  a  few  tens  of  m. 

The  footprint  of  a  single  pixel  is  shown  schematically  in  figure  2.  Prior  to  the  work  of 
Wollenweber  (1988),  LOWIllAN  6  calculated  the  radiance  along  the  atmospheric  path  from  the 
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camera  to  the  footprint.  This  atmospheric  path  is  labeled  “path  radiance'’  in  figure  2.  This 
calculated  radiance  was  taken  as  the  total  radiance  reaching  the  camera.  Infrared  radiance 
measurements  of  the  sea  and  sky  in  the  vicinity  of  the  ocean  horizon,  however,  often  resulted  in 
larger  radiances  than  those  predicted  by  LOW^^RAN  6,  pointing  out  the  need  for  additional 
radiance  mechanisms  in  such  a  setting. 

IWo  addib^.ial  radiance  contributions  introduced  by  Wollenweber  (1988)  can  be  understood 
as  follows.  The  pixel  footprint  will  contain  numerous  individual  wave  facets,  each  one  perhaps  a 
mm  or  less  in  extent  and  each  one  with  its  own  orientation,  which  constantly  changes  with  time 
as  the  wind  ruffles  the  open  sea  surface.  Within  the  footprint,  a  focet  with  a  given  instantaneous 
slope  will  reflect  a  portion  of  the  sky  dome  into  the  camera.  This  radiance  is  labeled  “sky  dome 
radiance”  in  figure  2.  Furthermore,  the  ocean  itself  radiates  like  a  blackbody  at  the  temperature 
of  the  water,  so  each  facet  will  emit  blackbody  radiation  into  the  sky  dome  above  it  and  some  of 
the  radiation  will  enter  the  camera.  This  radiance  is  labeled  “emissions”  in  figure  2.  These  two 
additional  sources  of  infrared  radiance  were  added  to  the  LOWTRAN  6  computer  code  by 
Wollenweber  in  1988.  Richter,  Hughes,  and  Paulson  (1989)  showed  that  these  maritime 
additions  to  LOWTRAN  6  substantially  increased  the  agreement  between  predicted  and 
observed  radiances  near  the  sea  horizon  in  the  8-  to  12-micron  spectral  region. 

At  this  point,  we  should  note  that  in  this  wavelength  band  (8  to  12  microns),  observations 
have  shown  that  glitter  is  not  a  significant  source  of  radiance.  Since  1989,  however,  rapid  strides 
in  detector  technology  in  the  3-  to  5-micron  spectral  region  have  intensified  interest  in  this 
spectral  band  and,  here,  sun  glint  is  known  to  make  a  major  contribution.  This  could  have  been 
expected  from  figure  1,  which  shows  two  trends  occurring  when  spectral  radiance  moves  from 
the  middle  of  the  8-  to  12-miCTon  band  (say  10  microns)  to  the  middle  of  the  3-  to  S-micron 
band  (say  4  microns):  First,  the  solar  spectral  radiance  rises  by  more  than  an  order  of  magnitude 
(the  actual  ratio  is  33)  in  going  from  10  to  4  pm  and,  second,  the  300-K  spectral  radiance  falls 
by  more  than  an  order  of  magnitude  (the  actual  ratio  is  15)  in  going  from  10  to  4  pm.  Hu^es 
(1992)  showed  that  even  with  the  two  Wollenweber  additions  for  reflected  sky  radiance  and  sea 
surface  emissions,  the  mode!  predictions  fall  short  of  the  sea  radiances  observed  in  the  presence 
of  sun  glint  (1992a). 

SUN  GLINT 

This  work  was  undertaken  following  the  suggestion  of  Hughes  that  sun  glint  be  included  in 
the  maritime  modifications  to  LOWTRAN  6  (1992a).  Some  indication  of  the  importance  of  solar 
reflections  to  infrared  shipboard  receivers  can  be  obtained  from  the  foil  .wing  statement: 
“Shipboard  defense  systems  employing  infrared  warning  receivers  have  experienced  saturation 
of  threat  detection  processing  circuits  due  to  the  signal  resulting  from  solar  reflection  (glints) 

from  the  rough  sea  surface - ”  The  analysis  shows  that  the  use  of  a  linear  polarizer  with  a 

shipboard  IR  warning  receiver  can  typically  be  expected  to  increase  the  number  of  daylight 
hours  that  the  sensor  can  be  operated  without  any  blanking,  and  to  narrow  the  azimuth  sector  for 

which  blanking  is  required  for  low  sun  angles - With  a  linear  polarizer  the  number  of  hours 

for  which  blanking  will  be  necessary  can  typically  be  reduced  to  6  to  8  hours  per  day  and  the 
azimuthal  extent  of  the  blanking  during  the  remaining  hours  reduced  to  ±  12*”’  (Beard,  1976, 
p.  5). 

The  seminal  paper  on  sun  glint,  or  glitter  as  it  is  sometimes  called,  is  by  Cox  and  Munk 
(1954).  This  paper  remains  definitive  to  this  day.  The  equations  and  results  of  Cox  and  Munk 
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together  with  the  calculations  available  within  the  LOWTRAN  6  code  will  be  the  basis  of  the 
work  presented  here. 

The  glint  has  been  included  according  to  the  following  approach.  For  a  perfectly  calm  sea, 
the  entire  pixel  footpirint  will  reflect  the  sun’s  rays  into  a  properly  oriented  camera.  An  eauatioc 
given  by  Cox  and  Munk  (1954)  can  be  used  to  relate  the  solar  spectral  irradiance  (W  cm*^  pm'^) 
falling  on  the  sea  to  the  spectral  radiance  (W  cm'^  str'^  pm*’)  leaving  the  sea.  When  multiplied 
by  the  spectral  extinction  transmittance  along  the  entire  path  from  footprint  to  camera,  this 
radiance  is  the  sun  glint  contribution  we  are  seeking. 

In  reality,  however,  only  a  certain  fraction  of  the  facets  will  be  correctly  positioned  within 
the  footprint  to  let  a  mirror-like  specular  reflection  of  the  sun’s  rays  into  the  camera.  Given  the 
position  of  the  camera  and  the  sun,  it  is  a  simple  matter  of  geometry  to  calculate  the  facet  slope 
required  for  a  specular  reflection.  Cox  and  Munk  (1954)  show  that  once  the  wind  velocity  is 
known,  the  fraction  of  facets  whose  slopes  lie  near  a  given  value  can  be  calculated.  This  spatial 
fraction,  which  can  also  be  regarded  as  a  temporal  probability  function,  is  then  applied  to  the 
radiance  described  above  in  order  to  provide  a  more  realistic  representation  of  the  sun  glint 
radiance  reaching  the  camera.  This  contribution  is  labeled  “solar  irradiance"  in  figure  2. 
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Figure  1  (overleaf).  Theoretical  spectral  radiance  predicted 
by  the  Planck  equation  over  p  for  ideal  blackbodies  of  6000  K 
(Ae  red  line,  ap}Hoximating  the  sun)  and  300  K  (the  blue  line, 
a|^>roximating  the  sea  and  sky).  The  blue  dashed  line  is  a 
nominal  value  for  the  radiance  of  the  sky  produced  by  diffusely 
scattered  siuilight  The  combined  effect  of  scattered  sunlight  and 
radiation  enait^  from  the  atmosphere  at  the  earth’s  temperature 
can  be  estimated  by  adding  the  two  blue  curves,  and  in  many 
situations  this  comes  close  to  die  actual  spectral  radiance  of  the 
daytime  sky  (Bell  et  al.,  1S>60). 
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Figure  2  (overleaf).  Schematic  representation  of  the  maritime 
contributions  to  LOWTRAN.  The  infrared  camera,  typically 
located  on  board  ship,  produces  a  picture  consisting  of  a 
sequence  of  picture  resolution  elements,  or  pixels.  The  pixel 
footprint,  typically  10  m  on  a  side,  consists  of  that  portion  of  the 
ocean  whose  radiance  contributes  to  the  formulation  of  a  single 
pixel.  Unmodified  LOWTRAN  calculates  the  radiance  along  the 
path  from  the  camera  to  the  pixel  footprint  (“path  radiance”).  The 
marine  modifications  take  into  account  the  following  additional 
contributions:  blackbody  radiation  from  the  sea  itself 
(“emissions”);  the  radiance  from  the  entire  sky  hemisphere 
reflected  in  the  surface  of  the  sea  (“sky  dome  radiance”);  and 
reflected  solar  irradiance,  or  sun  glint  (“solar  irradiance”). 
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SUN  GLINT  FORMULATION 


The  geometry  of  the  point  of  reflection  is  shown  in  figure  3,  The  coordinate  system,  its 
orientation,  and  the  meaning  of  all  the  angles  are  discussed  in  detail  in  appendix  A. 

If  the  instantaneous  slope  of  a  facet  in  the  upwind  direction  is  tana,  and  if  the  instantaneous 
slope  of  the  facet  in  the  crosswind  direction  is  tanP,  then  let  the  fraction  of  facets  whose  slopes 
lie  within  ±  (dtana)/2  of  tana  and  ±  (dtanP)/2  of  tan^  be  given  by  the  probability  function 


pftana,  tan/3,  W)dima  f/tan/5  , 


(1) 


where  it  has  been  anticipated  that  the  probability  will  depend  on  the  wind  velocity  W  (ms'^). 

If  the  spectral  solar  irradiance  ///"“'  (W  cm'^  falling  on  the  facet  can  be  found,  then 

the  spectral  radiance  (W  cm'^str^  leaving  the  facet  can  be  calculated  from  the 
following  formula  originally  given  by  Cox  and  Munk  (1954)  and  later  adapted  for  use  near  the 
horizon  by  Saunders  (1968): 


^  4cos‘*y  cos(d)  ^ 


Here,  R  is  the  reflectivity  of  the  facet,  Q  is  the  angle  of  incidence,  X  is  the  wavelength,  y  is  the 
facet  tilt,  6  is  the  zenith  angle  of  the  reflected  ray,’  and  S  is  the  Saunders  shadowing  factor.  The 

radiance  leaving  the  facet  must  be  multiplied  by  the  spectral  extinction  transmittance  Te(X,) 
along  the  entire  path  between  the  facet  and  the  camera  and  integrated  over  the  spectral  band  to 
obtain  the  final  result  for  the  contribution  of  the  glint  IV,  to  the  total  radiance  (W  cm'^  str'^) 
reaching  the  camera: 


^2 

AT  =  j  Nf°^^‘re(X)cL\  . 

A, 


(3) 


refraction  is  neglected,  angle  t  is  180®  minus  the  zenith  angle  of  the  camera, 
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Figure  3  (overleaf).  The  geometry  of  reflection  from  a  tilted 
wave  facet  on  the  ocean  surface.  The  plane  of  incidence  contains 
the  incident  ray,  the  facet  normal  (not  shown),  and  the  reflected 
ray.  The  x-y  plane  is  horizontal  and  the  x  axis  is  chosen  to  point 
upwind.  A  full  discussion  of  this  geometry,  including  the 
meaning  of  the  various  angles,  is  given  in  the  appendix. 
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DETAILS  OF  THE  FORMULATION 


The  following  will  explain  how  each  term  in  equations  (2)  and  (3)  was  obtained  and  provide 
examples  of  the  terms’  values  that  were  appropriate  to  the  situation  obtained  on  10  October  1991 
at  1450  PDT  (hereafter  called  “10  October”)  when  sun  glint  data  were  taken  at  Point  Loma  in 
San  Diego,  CA. 

Let  azimuths  q  be  measured  positive  clockwise  from  true  north  when  looking  down  and  let 
Qdetector  be  the  azimuth  of  the  direction  in  which  the  camera  is  pointing. 

Wind  velocity  W  is  taken  from  measurements  of  wind  speed  and  wind  direction.  The  wind 
direction  sets  the  orientation  of  the  x  axis  and  thereby  determines  the  relative  azimuths  of  the 
camera  (Qdetector  •  Qwind  =  0  +  Jt)  and  the  sun  (Qsource  -  Qwind  -  ^)»  with  respect  to  the  wind 
(see  appendix  A). 

The  facet  slopes,  tana  and  tan^,  and  the  cosine  of  the  angle  of  incidence,  Q,  are  calculated 
from  the  following  equations  whose  derivation  is  given  in  appendix  A: 


tana  =  (-  sin0cos0  -  sinpi  cosx)/(cos6  +  cos//) 
tan^  =  (+  sin0  sin0  +  sin//  sim)/(cos6  +  cos//)  (4) 

cosQ  -  [(1  +  sin0  cos0sin//  cosx  4-  sin0sin0sin//sinx  +  cos^  co^)/2]^/^ 


On  10  October,  the  wind  speed  was  westerly  (Qwind  =  90°)  at  4.6  m  s'^;  the  camera  was  centered 
on  the  horizon  (9  =  90°)  pointing  south  southwest  (Qwind  =  225°);  and  at  1450  PDT  the  solar 
zenith  angle  \i  was  50.47°  at  an  azimuth  Qsource  of  224.79°.  With  these  values,  equation  (4) 
gives  14.3°  for  a,  -14.3°  for  P,  and  70.2°  for  Q. 

It  is  easy  to  show  from  the  equations  given  by  Cox  and  Munk  (1954)  that  the  tilt  (y,  the 
angle  between  the  horizon  and  the  direction  of  steepest  ascent  within  the  facet)  is  related  to  the  x 
and  y  slopes  by 


tan^  =  tan^a  +  tan^;S  . 


(5) 


On  10  October,  the  tilt  was  19.8°. 

Cox  and  Munk  (1954)  determined  p{tana,  tanP,  W)  for  the  ocean  surface  by  means  of 
airborne  photographic  measurements.  They  found  that 
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(6) 


,  />  ^-1  1/  tan^a  tan^5\ 

pftana,  tai^,  W)  =  Qxajj^  *  exp  -  5I  +  ~^) 


ai  =»  0.000  +  3.16  •  lO-^W 
of  *  0.003  +  1.92  •  10-^iy 


where  Oy  and  are  the  standard  deviations  in  the  upwind  (x)  and  crosswind  (y)  directions, 
respectively.  Both  depend  on  the  wind  speed,  which  must  be  entered  into  these  expressions  in 
ms'^  A  wind  spee^  of  1  ms*^  is  typical  of  a  calm  sea;  a  wind  speed  of  10  ms'^  is  typical  of  a 
rough  sea. 

For  a  particular  wind  speed,  the  value  of  p  in  equation  (6)  will  oscillate  up  and  down  about 
an  average  value  q  given  by  equation  (6),  with  Ou^  and  Oc^  each  replaced  by  its  average  value 
a,v^  =  (Ou^  +  Oc2)/2: 


^^tany,  W)  =  (2nalvT^  exp 


02w  =  0.0015  +  2.54  •  lO-^iy 


(7) 


We  shall  call  q  the  nondirectional  wave  slope  probability. 

Of  all  the  factors  in  equation  (2),  p  is  the  most  rapidly  varying  function  of  geometry  and 
meteorology.  For  most  cases  of  interest,  p  falls  within  about  a  factor  of  2  of  q,  whereas  q  itself 
varies  between  100  and  0.01.  Therefore,  the  value  of  can  be  estimated  by  inspecting  the 
value  of  q',  in  the  3-  to  5-micron  region,  a  nondirectional  probability  of  0.1  or  more  is  required 
for  sun  glint  to  make  a  significant  contribution  (0.1  °C  or  more)  to  the  observed  radiance. 

Figure  4  shows  how  q  varies  with  solar  zenith  angle  (or,  equivalently,  solar  elevation)  for 
the  case  of  a  horizontal  view  (6  =  90°)  down  the  center  of  a  glint  pattern  toward  the  sun 
(Odetector  =  Osource)-  The  following  four  figures  (figures  5-8)  show  the  value  of  9  as  the  camera 
moves  off  the  center  of  the  glint  pattern.  Figures  5  to  8  refer  to  a  fixed  solar  zenith  angle  with 
the  camera  viewing  horizontally  at  various  wind  speeds.  Figure  8  (p  s  50°),  displaying  a  wind 
speed  of  5  ms  *  for  10  October,  shows  that  q  was  in  the  vicinity  of  0.1. 

The  shadowing  factor  was  introduced  by  Saunders  (1968)  to  account  for  the  fact  that  one 
facet  can  hide  another  when  the  sea  surface  is  viewed  at  very  shaUow  glancing  angles  on  the  far 
horizon.  Further  motivation  for  this  facto’’  was  the  unbounded  nature  of  equation  (2),  absent  5, 
near  the  horizon — that  is  for  cos(&)  in  the  denominator  of  equation  (2)  approaching  zero. 
“Without  the  inclusion  of  shadowing,  we  find  that  when  [cos(9)]  becomes  very  small,  the 
radiance  calculated  from  equation  (2)  becomes  unbounded;  by  including  5,  which  goes  to  zero 
with  [cos(0)]  we  ensure  that  the  radiance  calculated  from  equation  (2)  remains  finite”  (Saunders, 
1968,  p.  648).  Saunders’  shadowing  factor  is 
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S(e,  H9  =  2[1  +  erf(v)  +  (vjny^  expf- 

tan(>r/2  -  6) 

“  a 

<p-  *  2ff<r*  * 

=  0.003  +  5.12  •  10 , 


(8) 


where  erf{v)  is  the  error  function.  Figure  9  gives  the  value  of  the  ratio  5(9,  IV)/cos(9)  for  various 
wind  speeds.  Within  a  degree  of  the  horizon,  this  ratio  is  approximately  20  for  a  wind  speed  of 
5  ms'^  which  was  the  approximate  value  on  10  October. 


^  =  90* 

0  =  P 
■C.  I  A 


SO  85  80  75  70  65  60  55  50  45 
Solar  Zenith  Angle  n  (deg) 


Figure  4.  Nondirectional  wave  slope  probability  as  a  function  of 
solar  zenith  angle  for  a  horizontal  outlook  centered  on  a  glint 
pattern.  Probabilities  in  the  shaded  area  do  not  substantially 
contribute  to  the  radiance. 
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HORIZONTAL  VIEW  AT  SOLAR  ELEVATION  10 


0  5  10  15  20  25  30 

Camera  Azimuth  Referred  to  Sun  (deg) 


Figure  5.  Nondirectional  probability  as  a  function  of  camera 
azimuth  for  a  solar  elevation  of  10°.  The  camera  is  viewing  the 
horizon  and  the  pattern  is  symetrical  about  the  vertical  axis. 


HORIZONTAL  VIEW  AT  SOLAR  ELEVATION  20  ‘ 


100 


0  5  10  15  20  25  30 

Camera  Azimuth  Referred  to  Sun  (deg) 


Figure  6.  Nondirectional  wave  slope  probability  as  a  function  of 
camera  azimuth  for  a  solar  elevation  of  20°. 
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HORIZONTAL  VIEW  AT  SOLAR  ELEVATION  30 


0  5  10  15  20  25  30 

Camera  Azimuth  Referred  to  Sun  (deg) 


Figure  7.  Nondirectional  wave  slope  probability  as  a  function  of 
camera  azimuth  for  a  solar  elevation  of  30°. 


HORIZONTAL  VIEW  AT  SOLAR  ELEVATION  40* 


0  5  10  15  20  25  30 


Camera  Azimuth  Referred  to  Sun  (deg) 


Figure  8.  Nondirectional  wave  slope  probability  as  a  function  of 
camera  azimuth  for  a  solar  elevation  of  40°. 
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THE  SAUNDERS  SHADOWING  RATIO 


Figure  9.  The  Saunders  shadowing  factor  as  a  function  of  the  wind 
speed  and  zenith  angle  of  the  outgoing  (camera)  ray. 

The  reflectivity  of  the  sea  water  facet  R(Q,  X.)  is  given  by  Fresnel’s  equations  governing  the 
reflection  of  unpolarized  light  by  a  nonmagnetic  conducting  medium  with  a  complex  refractive 
index  a(X)  +  ip(X)  (Stratton,  1941;  pp.  505, 506;  eqs  74,  77): 

R(Q,X)  *  (Rp  +  R^/2 

^  ^  (q  -  cosD)^ 

^  (q  +  cosC)2  + 

^  _  (fa^  -  p^JcosQ  —  q)^  -f  (2aficosQ  -  p)^ 

^  ([o^  ”  P^JcosQ  +  q)2  +  (2aficosQ  +  p)^ 

+  ^  +  JiFTp)/!  (9) 


=  (  -  g  +  +  g^)/2 

g  =  a}  -  -  sio^Q 


h  =  2afi  . 


The  index  for  sea  water  between  3  and  28  was  taken  from  literature  (Hale  and  Querry, 
1973;  Querry  et  al.,  1977).  Figure  10  gives  the  value  of  the  reflectivity  of  sea  water  along  the 
center  of  a  glint  pattern  as  a  function  of  wavelength  for  a  variety  of  solar  elevations;  the 
reflectivity  was  l^tween  10%  and  15%  for  10  October. 

The  spectral  extinction  transmittance  and  the  spectral  solar  irradiance  (figure  11)  reaching 
the  facet  were  obtained  from  LOWTRAN. 


FORTRAN  CODING 

The  preceding  equations  were  introduced  into  the  LOWTRAN  6  FORTRAN  computer  code 
along  lines  established  by  Wollenweber  (1988).  For  those  familiar  with  LOWTRAN,  a  brief 
summary  of  Wollenweber ’s  approach  is  now  given. 

LOWTRAN  calculates  radiance  while  tracing  a  ray  backwards  from  the  detector  toward  the 
edge  of  the  atmosphere.  If  the  ray  remains  in  the  atmosphere,  then  no  alteration  is  made  to  the 
original  radiance  calculation  and  the  result  is  the  radiance  of  a  pixel  in  the  sky.  If,  on  the  other 
hand,  the  ray  hits  the  lower  boundary  of  the  atmosphere,  then  this  boundary  is  reinterpreted  as 
the  surface  of  the  sea,  the  calculation  is  interrupted,  and  the  radiance  computed  thus  far  is  saved 
as  the  path  radiance  of  figure  2.  During  the  interruption,  a  Wollenweber  subroutine  called  WAVE 
typically  issues  8  to  12  new  LOWTRAN -compatible  “geometry  cards”  representing  the  sky 
dome.  Each  new  card  has  height  0  and  a  2enith  angle  that  the  Cox-Munk  wave  slope  statistics 
have  determined  contributes  at  least  5%  to  the  overall  sky  radiance  reaching  the  camera  after 
reflection  in  the  sea  surface.  (The  geometry  of  this  situation  corresponds  to  the  geometry  of 
problem  one  in  appendix  A.)  The  calculation  is  then  resumed  to  derive  the  radiance  along  each 
of  the  new  sky  dome  paths.  All  sky  dome  radiances  are  reflected  in  the  sea  and  added  together  to 
arrive  at  the  contribution  labeled  “sky  dome”  in  figure  2. 

In  addition,  the  geometry  of  the  reflection  is  fully  determined  (problem  one  in  appendix  A 
again)  for  each  new  zenith  angle  that  has  been  introduced,  since  the  geometry  of  the  camera  ray 
is  fixed  beforehand.  Hence,  the  angle  of  incidence  is  also  known,  and  can  be  used  to  calculate 
the  emissivity  (1  minus  the  reflectivity')  of  the  facet,  thus  arriving  at  the  contribution  labeled 
‘  emissions”  in  figure  2. 


For  the  glint  calculation,  the  mode  is  temporarily  changed  from  “atmospheric  radiance” 
(lEMSCT  =  2)  to  “directly  transmitted  solar  inadiance”  (lEMSCT  =  3),  and  a  final  new 
“geometry  card”  is  issued  from  the  WAVE  compatible  with  the  new  solar  irradiance  mode. 

^j^facct  jg  iijgjj  obtained  upon  re.sumption  of  the  interrupted  LOWTRAN  calculation.  (The  solar 
position  has  already  been  introduced  by  means  of  a  new  card  [Card  2E]  in  the  LOWTRAN  input 
file.)  With  the  source  and  detector  positions  fixed,  the  geometry  of  the  reflection  is  once  again 
determined  (problem  two  in  appondix  A),  and  new  subroutines  carry  out  the  remainder  of  the 
calculations  required  by  equation  (2)  resulting  in  the  glint  contribution  labeled  “solar  irradiance” 
in  figure  2. 

If  the  original  ray  has  hit  the  lower  boundary  of  the  atmosphere,  then  the  radiance  contribu¬ 
tions  described  in  the  preceding  three  paragraphs  are  added  together  and  the  result  is  the  radiance 
of  a  pixel  in  the  sea. 
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HORIZONTAL  VIEW  TOWARD  SUN 


Wavelength  Otir.) 


Figure  10.  The  reflectivity  of  sea  water  as  a  function  of  optical 
wavelength  for  various  solar  zenith  angles.  The  view  is  horizontal 
along  the  center  of  the  glint  pattern. 

INCIDENT  IRRADIANCE.  10  OCTOBER  1991 


3  3.5  4  4  "  5 


Wavelength  (/im) 

Figure  11.  The  spectral  irr  diance  for  solar  radiation  incident  on 
the  ocean  wave  facet  at  1450  PDT  on  10  October  1991.  The 
shaded  area  is  6.889  x  10  W  cm 
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Due  to  the  general  organization  of  UDWTRAN,  equations  such  as  (3)  are  not  executed  for 
the  marine  modiflcations.  Equations  similar  to 

A]  A2 

^1 

are  executed  instead.  This  is  the  Achilles  heel  of  the  marine  modifications  to  LOWTRAN, 
because  there  is  always  the  possibility  that  an  unusual  set  of  atmospheric  conditions  or  ray 
parameters  will  cause  a  large  error  when  we  use  this  approximation.  However,  test  cases  show  an 
error  of  5%  to  10%  during  its  use  under  normal  conditions. 

COMPARISON  WITH  EXPERIMENT 

Measurements  of  the  near  horizon  radiance  over  the  open  ocean  were  made  in  the  infrared  at 
3  to  5  pm  and  8  to  12  pm  by  Hughes  with  a  calibrated  th  ermal  imaging  system  (AGA  THER¬ 
MOVISION  model  780)  (1992b).  The  system  has  a  2.9.'>“  field  of  view  and  an  instantaneous 
Held  of  view  of  0.87  mrad.  Simultaneous  measurements  of  meteorological  parameters  were 
acquired  on  the  surface  and  with  altitude  by  using  balloon-borne  radiosondes.  Measurements 
were  made  on  six  different  days,  during  wUch  temperature  inversions  occurred  at  different 
heights.  The  surface  winds  were  onshore  in  all  cases  (either  northwesterly  or  v  esterly)  at  low  to 
moderate  speeds,  the  horizon  was  visible  to  the  naked  eye,  and  the  cloud  coveiage  did  not 
exceed  scattered  conditions.  The  camera  was  set  to  a  fixed  position  at  a  given  altitude.  As  time 
passed,  solar  glints  (if  any)  were  brought  into  the  camera’s  Held  of  view.  Three  days  have  been 
selected  to  compare  the  calibrated  thermal  imaging  system  with  LOWTRAN  predictions.  These 
days  will  be  considered  in  order  of  the  inaeasing  importance  of  solar  effects  in  the  3-  to  5-pm 
band. 

14  MARCH  1991 

On  14  March  1991,  data  were  taken  at  0936  PST  with  the  camera  located  at  an  altitude  of 
23  m  and  directed  over  the  ocean  on  a  bearing  of  225°.  The  wind  was  westerly  at  6.2  ms'^ 

There  was  a  high  inversion  level  near  1900  m.  This  was  in  the  morning,  the  solar  azimuth  was 
129°,  the  solar  zenith  angle  was  49°,  and  the  effects  of  sun  glint  on  the  ocean  surface  were  not 
noticeable.  The  measured  3-  to  5-pm  radiances  are  shown  in  figure  12,  which  is  a  vertical 
radiance  slice  consisting  of  a  column  of  pixels  taken  from  the  full  iiriage.  The  box  in  the  upper 
right-hand  comer  of  the  figure  is  a  schematic  representation  of  the  a;:imuths  of  the  wmd  (W), 
sun  (S),  and  camera  (C)  at  the  time  the  data  were  taken.  The  small  vertical  arrow  indicates  true 
north.  In  the  main  part  of  the  figure,  the  region  above  90.15  represents  the  sky  and  the  region 
below  90.15  represents  the  sea.  On  this  day  there  was  almost  no  difference  between  sea  and  sky 
radiances. 

This  raises  the  question  of  how  the  horizon  is  assigned  in  the  infrared  image.  The  assignment 
is  IT  de  at  a  very  distinctive  cusp  in  the  8-  to  12-pm  data  near  the  visible  horizon,  then  retained 
for  the  3-  to  5-pm  band.  Furthermore,  in  the  8-  to  12-pm  band,  the  calculated  radiance  is  forced 


■:t(k)dX 


(10) 


to  match  the  measured  radiance  at  the  horizon  by  empirically  adjusting  the  continental  air  mass 
factor  (ICSTL  in  LOWTRAN)  and  the  visibility  (VIS  in  LOWTTIAN).  These  parameters  are 
then  retained  unchanged  for  calculations  in  the  3-  to  5-pm  band,  where  we  generally  discover 
that  measured  and  ralculated  radiances  obtained  by  these  procedures  do  NOT  agree.  This 
disagreement  will  be  much  more  evident  in  the  subsequent  figures. 

The  total  radiance  calculated  for  the  sky  and  sea  is  shown  after  adjustment  at  the  horizon 
according  to  the  procedure  of  the  previous  paragraph  by  the  solid  line  in  figure  12.  Although  a 
close  agreement  between  the  calculation  and  the  measurement  appears  at  first,  the  discrepancy  is 
actually  large,  on  the  order  of  l^C — a  value  that  should  be  compared  to  the  measurement 
accuracy  of  several  tenths  of  a  degree  C.  The  total  computed  radiance  has  been  broken  down  in 
figure  12  into  the  path,  sea,  and  sky  contributions.  There  is  no  glint  component,  which  is  not 
surprising  if  we  recall  the  conditions  under  which  these  data  were  taken.  Without  the  marine 
modifications  introduced  by  Wollenweber,  the  calculated  radiance  would  follow  the  solid  line  in 
the  sky  and  the  dashed  line  in  the  sea.  It  is  obvious  that  the  modifications  have  moved  the 
calculations  much  closer  to  the  measurements.  The  same  is  mie  for  all  the  following  cases. 


VERTICAL  RADIANCE  SLICE:  14  MARCH  1991 


RADIANCE  (mW  cm  *  sir  ’) 

Figure  12.  Observed  and  calculated  3-  to  S-pm  radiance  values  for 
0936  PST  on  14  March  1991.  The  horizontal  line  at  90.15®  divides 
the  data  into  those  from  the  sky  (above  the  line)  and  from  the  sea 
(below  the  line). 

26  JULY  1990 

The  cameras  were  located  at  a  33'm  altitude  and  pointed  due  west.  The  wind  was  out  of  the 
west  at  1.8  m  s'^  beneath  scattered  cumulus  clouds.  A  strong  temperature  inversion  near  400  m 
was  recorded  by  the  radiosonde  launched  at  1428  PDT.  At  1536  PDT,  the  first  data  (figure  13) 
were  taken,  the  solar  azimuth  and  zenith  angle  were  259.4®  and  37.9®  respectivly,  and  no  glint 
effects  were  observed  visually.  At  1613  PDT,  the  last  data  (figure  14)  were  taken,  the  solar 
azimuth  and  zenith  angle  were  265.6°  and  45.6®  respectively,  and  the  3-  to  5-|im  sensor  was 
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VERTICAL  RADIANCE  SLICE:  26  JULY  1990 
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Figure  13.  Calculated  and  measured  radiance  values  in  the  3-  to 
5-fxm  band  on  26  July  1990  before  sola:  glints  were  visible  in  the 
camera  field  of  view. 


VERTICAL  RADIANCE  SLICE:  26  JULY  1990 
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Figure  14.  Calculated  and  measured  radiance  values  in  the  3-  to 
5-(im  band  on  26  July  1990  when  strong  visible  glint  effects  were 
filling  the  camera  field  of  view.  There  is  no  Co.x-Muak 
contribution. 
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saturated  as  the  visually  observed  glint  filled  the  camera  field  of  view.  In  the  meantime,  the  8-  to 
12-fjim  sensor  did  not  saturate  and  remained  unchanged  between  the  first  and  last  data.  This 
agrees  with  expectations  that  solar  glint  effects  should  be  much  more  prominent  in  the  3-  to 
5-pm  band  than  in  the  8-  to  12*pm  band. 

What  is  surprising  is  that  the  Cox-Munk  model  shows  no  glint  contribution  for  any  of  these 
26  July  data!  In  fact,  at  1613,  the  calculated  glint  contribution  was,  at  most,  5  x  10'^  W  cm'^ 
str^  Apparently  not  enough  facets  were  present  with  the  slope  required  for  a  specular  reflection. 
Figure  8  shows  that  even  under  the  more  favorable  glint  conditions  of  a  50°  zenith  angle  and 
3  ms'^  wind  speed,  the  nondirectional  wave  slope  probability  is  in  the  vicinity  of  0.01,  which  is 
an  order  of  magnitude  too  low  for  an  observable  contribution  in  the  3-  to  5-pm  band. 

Without  any  glint  contribution,  the  calculated  radiance  values  are  identical  in  figures  13  and 
14,  because  they  are  based  on  the  same  atmospheric  data  and  viewing  geometry.  Only  the  solar 
position  has  changed.  Something  indeed  happened,  however,  because  a  large  change  occurred  in 
the  data  measured  between  1536  PDT  and  1613  PDT:  the  sky  radiance  rose  about  2°C  and  the 
sea  radiance  rose  about  5°C.  At  this  point,  we  must  conclude  that  either  the  Cox-Munk  model 
does  not  apply  or  that  there  is  some  other  solar  contribution  we  have  yet  to  take  into  account. 

10  OCTOBER  1991 

The  cameras  were  mounted  at  23  m  on  a  bearing  of  225°.  The  day  was  very  hazy  without 
clouds  The  wind  was  onshore  at  4.6  m  s*^  A  very  strong  temperature  inversion  began  at  the  low 
altitude  of  about  60  m  and  ended  at  428  m,  where  the  temperature  was  35  °C.  The  first  data  were 
taken  at  1153  PDT,  and  a  second  set  were  taken  at  1240  PDT,  when  the  sun’s  position  was  zenith 
angle  39.4°  and  azimuth  181°;  visible  sun  glint  effects  began  to  appear  at  1400  PDT,  and  the 
final  data  were  taken  at  1450  PDT  as  the  glint  pattern  was  exactly  centered  within  the  camera 
Held  (the  solar  position  at  1450  was  zenith  angle  50.5°,  azimuth  224.5°). 

The  measured  data  in  figure  15  show  a  slight  rise  in  the  sky  radiances  compared  to  the  sea 
radiances.  This  is  a  feature  that  Hughes  attributed  to  the  hot  atmosphere  close  to  the  sea  surface, 
which  is  due  to  the  strong  low  inversion  layer  (1992b).  The  calculated  and  measured  radiances 
match  in  the  sea  but  not  in  the  sky;  the  lack  of  a  match  in  the  sky  has  been  encountered  pre¬ 
viously. 

By  1450  PDT,  both  sea  and  sky  radiances  rose  substantially  and  the  sea  radiance  exceeded 
the  sl^  radiance.  (The  drop  in  measured  radiance  between  90.6°  and  91  °  occurred  in  the  kelp 
beds  about  1  km  off  the  Point  Loma  shoreline;  the  ocean  surface  in  these  kelp  beds  is  always 
smoother  than  the  surface  of  the  open  sea,  probably  because  of  the  effect  of  oil  released  by  the 
kelp)  (figure  16).  The  Cox-Munk  model  now  makes  a  contribution  to  the  radiance  comparable 
to  the  sea  emission  and  the  reflected  sky  radiance.  The  total  radiance  matched  the  shape  of  the 
measured  data  (apart  from  the  drop  in  the  kelp  beds),  but  a  large  offset  of  about  6°C  took  place 
in  both  sea  and  sky.  In  figure  17,  the  solid  line  labeled  “shifted”  was  derived  from  the  dotted  line 
labeled  “glint”  by  the  arbitrary  addition  of  0.37  mW  cm*^  str^  It  can  be  seen  that  the  shifted  fit 
is  extremely  good  outside  the  kelp  beds.  The  glint  is  necessary  for  this  fit;  the  dashed  line  in 
figure  17  shows  what  the  total  unshifted  radiance  would  be  without  glint.  It  is  interesting  to  note 
that  a  constant  offset  to  the  solid  line  of  figure  14  would  also  produce  a  calculation  running 
through  the  middle  of  the  data.  In  this  case,  the  offset  would  have  a  different  value  and  the  glint 
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VERTICAL  RADIANCE  SLICE:  10  OCTOBER  1991 
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Figure  15.  Calculated  and  measured  radiance  values  on  10  October 
1991  in  the  3*  to  5-^m  band  before  solar  glints  were  evident  in  the 
visible  spectral  region. 
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Figure  16.  Calculated  and  measured  radiance  values  on  10  October 
1991  in  the  3-  to  5-fim  band.  The  fall  off  in  the  measured  data  is 
due  to  the  Pt.  L>oma  kelp  beds. 
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VERTICAL  RADIANCE  SLICE:  10  OCTOBER  1991 
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Figure  17.  The  data  (dots)  and  total  radiance  (solid  line)  of 
figure  16  showing  the  effects  of  glint  removal  (dashed  line)  and  a 
constant  shift  (dotted  line). 


would  not  be  necessary  for  the  fit.  Judging  from  just  these  two  days  of  data,  an  additional 
mechanism  seems  to  exist  aeating  a  large  increase  in  both  sky  and  sea  radiance  at  just  the  time 
when  solar  glitter  (which  of  course  affects  the  sea  radiance  alone)  has  its  maximum  impact.  Such 
a  mechanism  might,  for  example,  be  multiple  solar  scattering,  but  it  would  have  to  enter  into 
both  sea  and  sky  with  equal  strength,  which  appears  to  indicate  that  it  must  be  independent  of 
path  length. 


CONCLUSION 

In  this  study,  wc  have  compared  3-  to  S-pm  radiance  data  measured  over  the  open  ocean 
(Hughes,  1992b)  to  the  radiance  predicted  by  the  LOWTRAN  6  computer  code  (Kneizys  et  al., 
1983),  which  was  modified  to  take  marine  effects  into  account  (Wollenweber,  1988).  The  marine 
effects  are  the  following: 

1.  Blackbody  emission  from  the  sea. 

2.  Radiance  from  the  sky  reflected  in  the  sea. 

3.  Solar  reflection  by  the  sea. 

The  data  were  taken  with  the  intent  of  revealing  radiance  contributions  due  to  the  sun  in  the 
3*  to  5-|im  band.  We  examined  five  sets  of  data  on  three  different  days  and  compared  the  data 
with  the  modified  LOWTRAN  6  predictions. 
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We  draw  the  following  conclusions  for  the  3-  to  S*(xm  band; 

1.  All  of  the  marine  modifications  make  an  important  improvement  in  the  fit  of  the 
mode!  to  the  data. 

2.  Significant  discrepancies  on  the  order  of  several  degrees  C  remain  between  the  model 
and  the  data. 

3.  Solar  glints  play  a  role  but  do  not  overwhelm  (or  even  dominate)  the  data. 

4.  Another  solar  effect  appears  to  play  as  strong  a  role  as  sun  glitter.  Its  origin  is 
unknown;  we  speculate  that  it  may  be  multiple  solar  scattering. 
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APPENDIX  A 


DERIVATION  OF  GEOMETRIC  EQUATIONS  GOVERNING 
THE  REFLECTION  OF  RADIATION  INTO  A  DETECTOR 
FROM  A  TILTED  OCEAN  WAVE  FACET 


THE  GEOMETRY 

The  geometry  at  the  point  of  reflection  is  shown  in  figure  3  in  the  main  text  (after  Cox  & 
Munk,  1954).  The  various  quantities  in  figure  3  as  well  as  others  that  will  be  required  are 
defined  by  the  following  concepts: 

The  Coordinate  System 

Choose  a  coordinate  system  whose  origin  is  the  point  of  reflection.  Let  the  x  axis  point 
upwind,  the  z  axis  toward  the  zenith,  and  the  y  axis  crosswind  such  that  a  right-handed  system  is 
lormed.  The  x-y  plane  is  tangent  to  the  earth  (horizontal)  at  the  point  of  reflection.  Let  azi¬ 
muthal  angles  be  measured  in  the  x-y  plane  (the  horizon)  and  consider  them  to  be  positive 
when  measured  clockwise  from  north  looking  down  along  the  -z  direction.  All  azimuths  point 
away  from  the  origin  except  the  azimuth  of  the  detector,  which  points  toward  the  origin  along 
the  line  of  detector  sight. 

The  Tilted  Facet 

Let  the  tilted  facet  pass  through  the  origin.  Let  tan  a  be  the  facet  slope  dz/dx  in  the  x-z 
plane,  and  let  tan  3  he  the  facet  slope  dzidy  in  the  y-z  plane.  Define  a  unit  vector  u„  normal  to 
the  facet  with  components  (n„,  Cn). 

The  Source 

The  incoming  ray  is  bent  due  to  the  effects  of  refraction.  At  the  point  of  reflection  define  p 
to  be  the  zenith  angle  of  the  ray  arriving  from  the  source  or,  in  other  words,  the  apparent  angle 
from  the  zenith  to  the  source  regarded  from  the  origin.  Azimuthal  angles  are  not  affected  by 
refraction.  Let  x  be  the  azimuth  of  the  source  with  respect  to  azimuth  of  the  wind: 

^  *  Qsource  ~  Qwind- 

Let  {%  be  a  unit  vector  (a,  ,  bj,  q)  along  the  incident  ray  pointing  from  source  to  origin. 

The  Detector 

The  outgoing  ray  is  also  bent  in  the  vertical  direction  due  to  the  effects  of  refraction.  Let  0 
be  the  zenith  angle  at  the  point  of  reflection  of  the  outgoing  ray,  or,  equivalently,  the  apparent 
zenith  angle  of  the  detector  as  observed  from  the  origin.  Let  0  be  the  azimuth  of  the  outgoing 
ray  with  respect  to  the  azimuth  of  the  wind  at  the  point  of  reflection: 

4*  =  Qdeiector  ~  Jt  ~  Qwind* 

Let  IV  be  a  unit  vector  (or,  br,  q)  along  the  reflected  ray  pointing  from  origin  to  detector. 
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The  Plane  of  Reflection 


The  plane  of  reflection  will  contain  Ui,  14,,  and  iv<  According  to  the  law  of  reflection, 


u,  -  U(  ^  2  co&Q  Uk  > 


(1) 


where  Q  is  the  angle  of  incidence  and  the  angle  of  reflection. 

PROBLEM  ONE 

A  major  contribution  to  the  signal  reaching  an  infrared  camera  viewing  the  ocean  is  the 
radiance  reflected  into  the  camera  from  various  portions  of  the  sky.  Each  ocean  facet  will  have 
its  own  tilt  at  any  given  moment  and  will  be  contributing  a  specular  reflection  from  a  unique  part 
of  the  sky  dome  into  the  camera.  In  this  situation,  we  want  to  know  what  part  of  the  sky  dome  is 
coming  into  play  for  a  given  facet  tilt  and  a  given  camera  orientation.  This  leads  us  to  our  flrst 
problem:  Assume  that  the  detector  orientation  (8, 4>)  is  fixed  with  respect  to  the  reflection  point. 
Choose  a  facet  slope  (a,  p).  Demand  a  specular  reflection  by  the  facet  from  some  portion  of  the 
sky  (the  source)  into  the  camera  (the  detector).  Problem  one  is  to  find  equations  for  the  angle  of 
incidence  (Q)  and  the  zenith  angle  of  the  source  (p)  in  terms  of  a,  P,  6,  and  (p. 

THE  VECTOR  COMPONENTS 

We  will  follow  the  approach  of  Cox  and  Munk  (1954)  and  first  find  the  components  of  the 
three  unit  vectors  that  appear  in  the  law  of  reflection. 

The  components  for  the  unit  normal  to  the  facet  are  found  by  noting  that  for  any  vector  r  = 

(x,  y,  z)  in  the  facet, 


tt*  •  r  =  fl;,*  +  +  c«z  =  0 


(2) 


with 


al +  b]i  + cl  =  1.  (3) 

By  evaluation  of  equation  (2)  in  the  x~z  and  y-z  planes,  the  ratios  a„/c„  and  b„/c„  are  determined, 
which  sets  the  direction  of  the  normal.  Then  c„  is  found  from  equation  (3),  which  sets  the  length 
of  to  one. 

The  components  of  the  incident  unit  vector  and  the  reflected  unit  vector  are  found  by 
resolving  them  first  into  components  normal  to  the  horizon  (parallel  to  the  z  axis)  and  within  the 
horizon,  then  resolving  the  horizontal  component  into  its  x  and  y  values. 

For  the  reflected  ray, 

Or  =  +  sin^cos^ 

br  ~  -  sin8  sin0  (4) 

Cr  =  +  cos^ 
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For  the  incident  ray. 


a,  =  —  sii^  cosx 
=  -4-  siiyi  sinx 
c,  =  -  CO^l 


(5) 


For  the  facet  normal, 

On  =  —  tana  wsy 

bn  -  -  tuifi  cosy  (6) 

c„  =  +  cosy 

where 

cos^  s  l/(tan^a  +  tan^^  +  1).  (7) 

As  shown  in  figure  1  in  the  main  text,  y,  the  tilt,  is  the  angle  formed  between  the  horizon  and  the 
direction  of  steepest  ascent  within  the  facet. 

THE  SOLUTION  TO  PROBLEM  ONE 

By  rewriting  the  law  of  reflection  as 

-  U/  ~  2cosQ  Un  -  Ur  (8) 

and  squaring,  we  find  that 

cosfl  =  Ur'Un  =  a^„  +  b,bn  +  C,Cn ,  (9) 

which  gives  an  expression  for  cos  C2  after  substitution  of  the  components  from  equations  (4)  and 
(6).  The  expression  for  cos  p  is  given  by  the  z  component  of  equation  (8)  after  substitution  of 
the  components.  The  results  are 

cosCi  ^  (-  sin0  cos/p  tana  +  sin9  sio0  taq^  +  cos9)(cosy) 

cosfi  *e  (-2sin6>  cos0  tana  +  2sinfl  sin^  taq^  -  [sec^  “  2]costf)(cosV)  (10) 

sccV  *  tan^a  +  tan^)9  +  1  =  tan^  +  1 . 

PROBLEM  TWO 

Another  major  contribution  to  the  camera  signal  will  be  sun  glint,  wherein  a  solar  ray  is 
reflected  into  the  camera  by  a  wave  facet.  In  this  situation,  we  know  the  source  and  detector 
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positions  but  need  to  know  the  wave  slope  required  for  a  specular  reflection.  This  brings  us  to 
the  next  problera:  Assume  that  the  detector  orientation  (6, 4>)  and  the  source  orientation  (p,  k) 
are  both  fixed  with  respect  to  the  point  of  reflection.  Demand  a  specular  reflection  by  the  facet 
from  the  source  into  the  detector.  Problem  two  is  to  find  the  angle  of  incidence  (Q)  and  the  facet 
slopes  (tana,  tan^)  in  terms  of  6,  <)>,  p,  and  h. 

THE  SOLUTION  TO  PROBLEM  TWO 

To  find  the  angle  of  incidence,  square  equation  (8)  to  obtain 

2cOS^f?  =  i  -  14,  •  u,  (11) 


and  substitute  component  values  for  reflected  and  incident  rays  from  equations  (4)  and  (5), 
respectively.  The  facet  slopes  are  found  by  evaluating  the  components  of  Un  in  equation  (1)  by 
using  equations  (4)  and  (5)  once  again.  The  result  is: 

cosf?  -  [(1  +  sinflcos^siqucosx  +  sindsin^siiiusinx  +  cos6co^)/2]*/^ 

tana  =  (-  sindcos^  -  siiV4COsx)/(cosd  +  cos;u)  (12) 

taq^  =  (+  sindsin^  +  siiV4sinx)/(cos^  +  cos^) 

tanV  *  tan^a  +  tan^y?. 

For  viewing  toward  the  sun,  <{>  *  x  -  n  and  the  above  equations  reduce  to 


a 

tana 

tai^ 

y 


pi  6 
2 

+  cos^tany 
-  sin^tany 


(13) 


Some  useful  additional  relationships  are 

cosf?  =  [(1  +  sin0siiv4cos(0  -  x)  +  cos^cosiu)/2]*/^ 
cos(2f2)  =  sin^siip4cos(0  ~  x)  +  cosdcospi  (14) 

cosy  =  (cosO  +  cospi)/2cosQ . 


SYMMETRIES 

The  signs  of  the  components,  hence  the  signs  of  the  terms  in  equations  (10)  and  (12),  depend 
on  what  conventions  are  employed  when  the  problem  is  setup.  TWo  changes  in  the  initial  setup 
that  might  be  considered  are  the  following: 

1.  Change  the  sense  of  azimuth  to  one  that  is  more  natural  for  the  x-y-z  coordinate  sys¬ 
tem;  namely,  measure  azimuthal  angles  positive  going  counterclockwise  looking 
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down  on  the  bori2on  along  the  -2  direction.  If  this  convention  is  adopted,  all  equa¬ 
tions  for  the  components  remain  the  same  except  that  the  sign  of  the  >^omponent  is 
reversed  for  the  incident  and  reflected  rays: 

br^-br 

^  A, 

t,  —  -  b, 

2.  Change  the  definition  of  ()>  to  Qdetector  -  Qwind  so  that  it  gives  the  angle  between  the 
direction  in  which  the  camera  is  pointed  and  the  direction  of  the  wind.  If  this  con¬ 
vention  is  adopted,  then  the  equations  for  all  components  remain  the  same  except  that 
the  signs  for  the  x  and  y  components  of  the  reflected  ray  are  reversed: 

a,-*  -  Or 


br^-br 

3.  If  both  changes  are  adopted  at  die  same  time,  then 

Or-*  -  Or 

bi  -*  -  f>. 

and  the  solutions  to  problems  one  and  two  become  respectively 

cosf3  s:  Csindcos^taDa  +  sin^io^tai^  -h  co^Xcosy) 

com  =  (2sindcos^tana  +  2sin^siD^taiv9[secV  "  2]cos0)(cos^) 

seeV  =  tan^G  +  tan^^  +  1 


and 


cosf2  =  ([-sin^cos^sinucosx  -  sin^sin^siqGsinx  +  cos^cos;m]/2)'/^ 
tana  =  ("+  sin0sin0  -  siqucosx>/fsin5  +  co^^ 
taivS  =  (+  sindsin0  -  sity/sinx^/fcos^  -h 


(16) 


(17) 


(18) 


(19) 


COMPARISON  WITH  FORTRAN  CODE 

Problem  one  has  previously  been  considered  by  Wollenweber  (1988)  and  introduced  into 
LOWTRAN  6  as  a  set  of  subroutines  that  account  for  the  reflection  of  infrared  sky  radiance  by 
ocean  wave  facets  into  an  infrared  camera.  The  equations  for  cosQ  and  cosp  occur  just  after  line 
100  of  the  subroutine  GEOPRO.  The  following  table  associates  the  terms  in  the  FORTRAN 
code  with  the  symbols  used  here: 


1  AMEN 

SX 

SY 

CAL 

PSI 

COW 

OME 

PHI 

1  This  work 

f  scc^Y 

tana 

tan^ 

cos£2 

COS|X 

Q 

Ttn-e 

Furthermore,  the  code  apparently  adopts  both  of  the  changes  mentioned  in  the  previous 
section  since  CAL  and  COW  both  have  positive  first  terms. 
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APPENDIX  B 


LIST  OF  SYMBOLS 
Braces  {}  indicate  subscripts 


a 


9 


LATIN 

a 

{av} 

b 

c 

{C} 

C 

{e} 

H 

{i} 

{n} 

N 

N>. 

P 

q 

R 

{r} 

{s} 

S 

u 

{^} 

W 


X  component  of  a  unit  vector  u 
Average 

y  component  of  a  unit  vector  u 
z  component  of  a  unit  vector  u 
Camera;  also  Crosswind 
Degree  Centigrade 
Extinction 
Irradiance  (W  cm*^ 

Spectral  irradiance  (W  cm-^  (rm*^) 

Incident 

Normal 

Radiance  (W  cm-^  sir'^) 

Spectral  radiance  (W  cm'  -  str^  pm'^) 

Cox-Munk  directional  wave  slope  probability 

Cox-Munk  nondirectional  wave  slope  probability 

Reflectivity 

Reflected 

Sun 

Saunders  shadowing  factor 

Unit  vector 

Upwind 

Wind  speed  (ras*^) 
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{w} 


Wind 


X  Upwind  direction  in  the  horizon  at  the  point  of  reflection 

y  Crosswind  direction  in  the  horizon  at  the  point  of  reflection 

z  Direction  of  the  zenith  at  the  point  of  reflection 

GREEK 

a  Angular  slope  of  wave  facet  in  x  direction;  also 

Real  part  of  the  optical  index  of  refraction 

P  Angular  slope  of  wave  facet  in  y  direction;  also 

imaginary  part  of  the  optical  index  of  refraction 

Y  Tilt,  the  angular  slope  in  the  direction  of  steepest  facet  ascent  measured  with 
respect  to  the  horizon 

6  Zenith  angle  of  reflected  ray  at  the  point  of  reflection 

V  Normalized  elevation  of  reflected  ray  at  point  of  refleaion 

X  Azimuth  of  sun  with  respect  to  azimuth  of  wind 

X  Wavelength  (fim) 

p  Zenith  angle  of  sun  at  point  of  reflection 

Q  Azimuth,  measured  positive  counterclockwise  from  true  north  looking  down  on 

the  horizon 

^  Root  mean  square  deviation 

t  Transmittance 

Q  Angle  of  incidence 
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13  ABSTRACT  (Mnrnwn  MO  wordi) 
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predict'd  and  observed  radiwces  differed,  however,  oy  a  large  constant  with  an  equiv^ent  blackbody  temperature  of 
approA..uately  5°C.  Personnel  speculated  that  the  origin  of  this  additional  radiance  may  be  multiple  scattering  of  solar  radi- 

ation. 
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